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Abstract 

We discuss the (g ip p )d self-interacting scalar field theory, in the strong-coupling regime. We 
assume the presence of macroscopic boundaries confining the field in a hypercube of side L. We 
also consider that the system is in thermal equilibrium at temperature For spatially bounded 
free fields, the Bekenstein bound states that the specific entropy satisfies the inequality Ji < 2%R, 
where R stands for the radius of the smallest sphere that circumscribes the system. Employing the 
strong-coupling perturbative expansion, we obtain the renormalized mean energy E and entropy 

_ 2 

5" for the system up to the order (g ) p , presenting an analytical proof that the specific entropy 
also satisfies in some situations a quantum bound. Defining as the renormalized zero-point 
energy for the free theory per unit length, the dimensionless quantity £ = f and hi(d) and ^(gQ 
as positive analytic functions of d, for the case of high temperature, we get that the specific 
entropy satisfies ^ < 2nR £. When considering the low temperature behavior of the specific 

entropy, we have | < 2nR ^£ 1 ~ d . Therefore the sign of the renormalized zero-point energy can 

£ d 

invalidate this quantum bound. If the renormalized zero point-energy is a positive quantity, at 
intermediate temperatures and in the low temperature limit, there is a quantum bound. 
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1 Introduction 



There have been a lot of activities discussing classical and quantum fields in the presence of 
macroscopic boundaries. These subjects raise many interesting questions, since boundaries intro- 
duce a characteristic size in the theory. For example, in the field-theoretical description of critical 
phenomena, the confinement of critical fluctuations of an order parameter is able to generate long- 
range forces between the surfaces of a film. This is known as statistical mechanical Casimir effect 
P] [2] [3] [3]. These long-range forces in statistical mechanical systems are characterized by the 
excess free energy due to the finite-size contributions to the free energy of the system. It should be 
noted that the statistical mechanical Casimir effect is still waiting for a satisfactory experimental 
verification. On the other hand, the electromagnetic Casimir effect [5], where neutral and perfectly 
conducting parallel plates in vacuum attract each other, has been tested experimentally with high 
accuracy. The introduction of a pair of conducting plates into the vacuum of the electromagnetic 
field alters the zero-point fluctuations of the field and thereby produces an attraction between 
the plates [6] [7] [8] [9]. A still open question is how the sign of the Casimir force depends on 
the topology, dimensionality of the spacetime, the shape of bounding geometry or others physical 
properties of the system [10] [11] [12] [13]. We should emphasize that the problem of the sign 
of the renormalized zero-point energy of free fields described by Gaussian functional integrals is 
crucial for the subject that we are interested to investigate in this paper. 

Another basic question that has been discussed in this scenario, when quantum fields interact 
with boundaries, is about the issue that these systems may be subjected to certain fundamental 
bounds. One of these proposed bounds relates the entropy S and the energy E of the quantum 
system, respectively, with the size of the boundaries that confine the fields. This is known as 
the Bekenstein bound which is given by S < 2ir E R/hc, where R stands for the radius of the 
smallest sphere that circumscribes the system [2] [15] [16] [TF] [18]. Such bound was originally 
motivated by considerations of gravitational entropy, a consistency condition between black hole 
thermodynamics and ordinary statistical physics that could guarantee that the generalized second 
law of thermodynamics is respected, which states that the sum of the black-hole entropy and the 
entropy of the matter outside the black-hole does not decreases. For example, in a Schwarschild 
black-hole in a four-dimensional spacetime, the Bekenstein entropy, which is proportional to the 
area of the spherical symmetric system, exactly saturates the bound. When gravity is negligible, 
the bound must be valid for a variety of systems. 

Although analytical proofs of this quantum bound on specific entropy for free fields has been 
proposed in the literature, many authors in the past criticized the bound [19] [20] [21] [22] [23]. 
Deutsch claims that the quantum bound is inapplicable as it stands to non-gravitating systems, 
since an absolute value of energy cannot be observed, and also that for sufficient low temperatures, 
a generic system in thermal equilibrium also violates the entropy bound. Unruh pointed out that 
for system with zero modes, the specific entropy cannot satisfy any bound. Many of these criticisms 
were answered by Bekenstein and collaborators. The problem of the low temperature systems was 
answered in Ref. [23] and the problem of systems with zero modes was answered in the Ref. |25j . 
An strong argument used in one of these examples is based in the fact that the renormalized zero- 
point energy of some free quantum field could be negative. Some authors claim that, if we take 
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into account the boundaries responsible for the Casimir energy, it is possible to compensate their 
negative energy, yielding a positive total energy which respects the Bekenstein bound, although 
this is far from a simple problem 

In fact this last point that we pointed out, i.e., the objection raised against the violation of the 
Bekenstein bound for free fields, must be analyzed more carefully. The problem now confronting 
us is to prove that, although Casimir energy can be negative for some physical situation, the 
sum of the energy of the boundary and the Casimir energy will be positive. In other words, the 
contribution from the boundary would make the total energy of the system always positive for any 
configuration of the hyperplane when the Casimir energy is negative. How close the boundaries 
must be in such a way that the positive contribution from the rest mass of the boundary is always 
smaller than the modulus of the Casimir energy? Our conclusions from the above arguments is 
that, for example, in the case of the electromagnetic field, it is essential to construct a microscopic 
model where effects of dispersion and absorbtion must be taken into consideration. Many authors 
discussed the problem of quantization of the electromagnetic field in dispersive and absorptive 
linear dielectrics [27] [28] [29] [30] . The essential question that confronts us is the positivity of the 
total energy of any quantum system defined in a compact domain in any situation. The validity 
of the Bekenstein bound for configurations with negative Casimir energy depends on the answer 
for this last question. Our intention in this paper is not to study such deep and difficult question, 
introducing physically realistic boundary conditions, but only to discuss the situation of idealized 
mirror boundaries. 

We may observe that another quite important situation has not been discussed systematically 
in the literature, at least as far we known. A step that remains to be derived is the validity of the 
bound for the case of interacting fields, which are described by non-Gaussian functional integrals, 
at least up to some order of the perturbation theory. Nonlinear interactions can change dramati- 
cally the energy spectrum of the system and this might lead to the overthrow of the bound [31] |32j . 
The difficulties that appear in the implementation of this program in the presence of macroscopic 
boundaries are well known. In systems where the translational invariance is lacking it is much 
harder to compute the Feynman diagrams then in the unbounded space. Nevertheless a regulariza- 
tion and renormalization procedure can in principle be carry out in any order of the perturbative 
expansion [33]. See for example refs. [34J [35J [36] [37], where the perturbative renormalization 
were presented in first and second order of the loop expansion in the A<£> 4 self-interacting scalar 
field theory. We would like to stress that the renormalization program is implement in a different 
way from unbounded or translational invariant systems because surface divergences appear. 

The aim of this paper is to show for a given self-interaction field theory in which situations the 
specific entropy satisfies a quantum bound. For the answer of this important question, there are 
two different routes. The first one is to use the weak-coupling perturbative expansion. However, 
as we discussed, in the renormalization procedure appears surface divergences that forces one to 
introduce non-local counterterms. It is unclear for us how this affect the physical relevance of the 
results that can be obtained. There are some kind of problems for which the mean energy and the 
canonical entropy a system can easily be found for quantum fields defined in a simple connected 
bounded region. We can show that using the strong-coupling expansion one can evaluate the 
mean energy and the canonical entropy of the system in a regime in which quantum fluctuations 
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dominate. 

Therefore we study the (go (p p )d self-interacting scalar field theory in the strong-coupling 
regime. We assume the presence of macroscopic boundaries that confine the field in a hyper- 
cube of side L and also that the system is in thermal equilibrium with a reservoir. We present an 

_ 2 

analytic proof that, up to the order (go) p , the specific entropy satisfies in some situations a quan- 
tum bound. Defining ey as the renormalized zero-point energy for the free theory per unit length, 
£ = t and h\(d) and /^(cQ as positive analytic functions of d, for the case of high temperature, 
we get that the specific entropy satisfies the inequality J? < 2tiR £. When considering the low 
temperature behavior of the specific entropy, we have | < 27iR^^- d . We are establishing; a 

£ d 

bound in the strong-coupled system in the following cases: in the high temperature limit and if 
the renormalized zero point-energy is a positive quantity, at intermediate temperatures and also 
in the low temperature limit. 

In the weak-coupling perturbative expansion, the information about the boundaries can be 
implemented over the free two-point Schwinger function Go(m ;x — y) of the system. In the 
strong-coupling perturbative expansion, we have to deal with the problem of how the boundary 
conditions can be imposed. Let us briefly discuss the strong-coupling expansion in Euclidean 
field theory at zero temperature. The basic idea of the approach is the following: in a formal 
representation for the generating functional of complete Schwinger functions of the theory Z(V, h), 
we treat the Gaussian part of the action as a perturbation with respect to the remaining terms 
of the functional integral, i.e., in the case for the {go<p p )d theory, the local self-interacting part, 
in the functional integral. In the generating functional of complete Schwinger functions, V is the 
volume of the Euclidean space where the fields are defined and h(x) is an external source. We 
are developing our perturbative expansion around the independent-value generating functional 
Qo(h) [38] [39] [ID]. In the zero-order approximation, different points of the Euclidean space are 
decoupled since the gradient terms are dropped [H] [12] [13] [H]. 

The fundamental problem of the strong-coupling expansion is how to give meaning to the 
independent-value generating functional and to this representation for the Schwinger functional. 
One attempt is to replace the Euclidean space by a lattice made by hypercubes. A naive use of 
the continuum limit of the lattice regularization, where one simply makes use of the central limit 
theorem for the independent-value generating functional, leads to a Gaussian theory. A solution 
to this problem was presented by Klauder a long time ago [38] [10] [11] • The modification which 
allows us to avoid this limitation is a change in the usual definition of the measure in the functional 
integral, which possesses local translational invariance, by another one which is non-translational 
invariant. 

Let us remark that, in the strong-coupling regime, assuming that the source is constant, we 
can perform the perturbative expansion around a independent- value generating function, up to 

_ 2 

the order (g ) p, and it is possible to split \iaZ(V, h) in two contributions: one that contains only 
the independent- value generating function and other that contains the spectral zeta-function. 
Therefore, in order to obtain the thermodynamic quantities, one must proceed in two stages. 
First, one gives a operational meaning to the independent- value generating function; then, one 
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consistently implements the boundary conditions in the strong-coupling regime. Since we are 
working in first order of perturbation theory, to implement boundary conditions, we use the 
spectral zeta-function method jlS] [46] . For a complete review about the subject see for example 
refs. |47J [48J. Quite recently a very simple application of this formalism was presented [41?] . where 
it was considered an anharmonic oscillator in thermal equilibrium with a reservoir at temperature 
f3^ 1 . Using the strong-coupling expansion, it was found the mean energy in the regime A ^> u, 
up to the order where A and u are the coupling constant and the frequency of the oscillator 
respectively. 

The organization of the paper is as follows: In section II we discuss the strong-coupling expan- 
sion for the (go V ?p )d theory. In section III we discuss he free energy and the spectral zeta-function 
of the system. In section IV we show that it is possible to obtain in some situations a quantum 
bound in the considered model. Finally, section V contains our conclusions. In the appendix 
A we present the Klauder's result, as the formal definition of the independent-value generating 
functional derived for scalar fields in a (^-dimensional Euclidean space. In the appendix B we prove 
that the spectral zeta-function Cd(s) evaluated in the extended complex plane at s = vanishes. 
To simplify the calculations we assume the units to be such that h = c = kg = I. 



2 The strong-coupling perturbative expansion for scalar 



Let us consider a neutral scalar field with a (g ip p ) self-interaction, defined in a (^-dimensional 
Minkowski spacetime. The vacuum persistence functional is the generating functional of all vac- 
uum expectation value of time-ordered products of the theory. The Euclidean field theory can 
be obtained by analytic continuation to imaginary time allowed by the positive energy condition 
for the relativistic field theory. In the Euclidean field theory, we have the Euclidean counterpart 
for the vacuum persistence functional, that is, the generating functional of complete Schwinger 
functions. In a d- dimensional Euclidean space, the self-interaction contribution to the action is 
given by 



The basic idea of the strong-coupling expansion at zero temperature is to treat the Gaussian 
part of the action as a perturbation with respect to the remaining terms of the action in the 
functional integral. Let us assume a compact Euclidean space with or without a boundary, where 
the volume of the Euclidean space is V. Let us suppose that there exists an elliptic, semi- 
positive, and self-adjoint differential operator O acting on scalar functions on the Euclidean space. 
The usual example is O = (—A + m^), where A is the (^-dimensional Laplacian. The kernel 
K(m ; x,y) = K(m ; x — y) is defined by 



(g Q (f p ) d theory 




(i) 




(2) 
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Using the fact that the functional integral which defines Z(V, h) is invariant with respect to the 
choice of the quadratic part, let us consider a modification of the strong-coupling expansion. We 
split the quadratic part in the functional integral which is proportional to the mass squared in two 
parts; one in the derivative terms of the action, and the other in the independent value generating 
functional. The Schwinger functional can be defined by a new formal expression for the functional 
integral given by 

Z(V,h)=exp^Jd d x J d d y^K(m ,a;x-y)-^j Q (a,h) } (3) 

where Qo(a, h), the new independent value functional integral, is given by 

Q (a, h)=M J [dip] exp (j ' d d x(- l -o m\ <p 2 (x) - ^ if{x) + h(x) V (x)) J , (4) 

and the modified kernel K(m , a;x — y) that appears in Eq. Q, is defined by 

K(m , <r;x-y)= (-A + (1 - a)m 2 ) 5 d {x - y), (5) 

where a is a complex parameter defined in the region < Re (a) < 1. 

The factor M is a normalization that can be found using that Qo( cr ) ^)U=o = 1- Observe 
that the non-derivative terms which are non-Gaussian in the original action do appear in the 
functional integral that defines Qo(a,h). At this point it is convenient to consider h(x) to be 
complex. Consequently h(x) = Re(/t) + ilm(h). In the paper we are concerned with the case 
Re(h) = 0. 

Since we are assuming a spatially bounded system in equilibrium with a thermal reservoir at 
temperature the strong-coupling expansion can be used to compute the partition function 
defined by Z(f3,Q,h)\ h=o, where h is a external source and we are defining the volume of the 
(d— 1) manifold as V^_i = Q. From the partition function we define the free energy of the system, 
given by F(f3,Q) = — ^ In Z(j3,Q,h)\h=o- This quantity can be used to derive the mean energy 
E(f3,Q), defined as 

E((3,n) = --^\nZ(f3,n,h)\ h=0 , (6) 

and the canonical entropy S(/3, Q) of the system in equilibrium with a reservoir with a finite size 
given by 

S((3,n) = (l-(3-^j\nZ(P,n,h)\ h=0 . (7) 
In the next section we will show that in a particular situation it is possible, up to the order 

_ 2 

(go) p to split In Z(P,n, h) in two parts: the first one that contains only the independent-value 
generating function and the second one that has the information on the boundary condition and 
it is given by derivative of the spectral zeta-function defined in the extended complex plane in 

8 = . 
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3 The independent- value generating function and the spec- 
tral zeta-function 



We are interested in global quantities. For simplicity we are assuming that the external source 
h(x) is constant. In this situation we call Z(V, h) as a generating function. Since we are in- 
troducing boundaries in the domain where the field is defined, the spectrum of the operator 
D = (—A + (1 — <j)rriQ ) has a denumerable contribution, and an analytic regularization proce- 
dure can be used to control the divergences of the theory. In order to impose boundary conditions 
the functional integral must be taken over functions restricted to the geometric configurations. At 

_ 2 

zero temperature, in the leading-order approximation (up to the to the order (g ) p) we can write 
the logarithm of the generating function as 

k m h) = <wU w h) (-! + \l &(s)l -) • (8) 

where a is a infinite constant and Cd( s ) is t ne spectral zeta-function associated with the elliptic 
operator D. 

Let us consider now the situation in which the system is finite along each one of the spatial 
dimensions, i.e., Xi G [0,L], i = 1,2, ...,d— 1. For the Euclidean time we assume periodic 
boundary conditions (Kubo-Martin-Schwinger KMS [50] [51J conditions) and for the Euclidean 
spatial dimensions we assume Dirichlet boundary conditions. We call this latter situation "hard" 
boundaries. See for example the Ref. [52] • For different kinds of confining boundaries see [53] 
[54] . The choice of the hard boundary provides an easy solution to the eigenvalue problem, so 
that explicit and complete calculation using the spectral-zeta function can be performed without 
difficulty. 

It follows that the operator D has the spectrum given by X ni ,...,n d where 



Tii7l\ 2 /nd_l7T\ 2 /2nrf7T 



(9) 



7ii,ri2, ... , Tid-i are natural numbers different from zero, since we are choosing Dirichlet boundary 
conditions and n,i are integer numbers. The spectral zeta-function associated with the operator 
D in this situation reads 



oo 



Cd(s)= E a,; „„. (10) 

ni,...,n d 



where s is a complex parameter, and the prime sign means that the term ri\ = 0, ri2 = 0, .., nd-i = 
must be excluded. The series above converges for Re s > \ and its analytic continuation defines 
a meromorphic function of s, analytic at s — 0. Since we should have to introduce an arbitrary 
parameter /i with dimension of a mass to implement the analytic procedure with dimensionless 
quantities, we have scaling properties. 
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Using n as a general index instead of rix, n^, the scaling properties follows from the fact 
that 

oo oo 

CM") = £( \ n y s = ^ E K s = » 2s Cn(s) . (ii) 

n n 

Therefore we have 

2~^^ d( ^' s=0 = ^Cd(s)| s=0 + ^ln/i 2 CD(s)| s=0 . (12) 

Before continuing, we would like to discuss two points. The first one is the fact that for different 
boundary condition, as, for example, Neumann boundary conditions in all the hyperplanes or 
periodic boundary conditions in all the spatial directions, the presence of the zero-mode can make 
the calculations more involved. Studying a two-dimensional spacetime T x S 1 and also T x S 3 , 
Dowker [55] discussed how the spatial zero mode contributes to the free energy. The zero-mode 
problem was also discussed by Dowker and Kirsten [56J. Elizalde and Tort [57] discussed the 
contribution to the thermal energy coming from the spatial zero mode in a system defined in a 
manifold with non-trivial topology; being more specific these authors studied a massive scalar 
field in S 1 x S 3 geometry. Kirsten and Elizalde discussed the Casimir energy of a massive scalar 
field in a general (2 + 1) - dimensional toroidal spacetime. For the massless case they excluded the 
n i — n 2 — mode [58]. They discussed also in an ultrastatic <i-dimensional spacetime how the 
ambiguity of the Casimir energy is related to the B±. The entropy bound related to various fields 

in the R x S 3 geometry was also discussed by many authors. See for example the Refs. [59] [60] . 
It is important to remark that this zero mode problem does not appear in the calculations that we 
are presenting, since we are choosing Dirichlet boundary conditions in all hyperplanes, excluding 
the possibility of the spatial zero mode. The generalized Bekenstein bound for systems with the 
spatial zero mode must be taken in consideration to extend the results of the paper. 

The second point is that it is possible to show that there is no scaling in the situation that 
we are interested in. The spectral zeta-function is related to the heat-kernel or diffusion operator 
via a Mellin transform. The trace of the diffusion operator is the integral of the diagonal part 
of the heat-kernel over the manifold. It is possible to perform an asymptotic expansion for the 
heat-kernel and this asymptotic expansion shows that the spectral zeta-function is a meromorphic 
function of the complex variable s possessing simple poles where the residues of the poles depends 
on the B n coefficients which depends on the Seeley-DeWitt coefficients, the second fundamental 
form on the boundary and the induced geometry on the boundary. See for example the Ref. [61] 
|62j . It is possible to show that the polar structure of the analytic extension of the spectral zeta 
function in a compact manifold with boundary is given by 



Cd(s) 



(47r)ir(s 



OO 75 



n=0 n 2 + S 



+ 92(s) 



(13) 



for n integer or odd-half integer, where g2{s) is an analytic function in C. As was stressed by 
Blau et al [63], in a four dimensional flat spacetime with massless particles and thin boundaries 
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the geometric coefficient B 2 vanishes. The spectral zeta function has poles at s = 1, s = 2, and 
so on. There would be a pole at s = 0, but is cancelled out by the pole in the gamma function. 
Therefore Cd( s )|s=o = Bd. For the case of hypercube with Dirichlet boundary conditions we prove 
in the appendix B that the spectral zeta-function in s = is zero, consequently Bd = and there 

is no scaling in the theory. We shall come back to this point in the conclusions. 

Let us study in Eq. ([S]) the contribution arising from the spectral zeta-function which takes 
into account the geometric constraints upon the scalar field. Using the spectrum of the D operator 
given by Eq. §§§ and the definition of the spectral zeta-function given by Eq. ffTUj) . we get that 
the derivative of the spectral zeta-function in s = yields 

5 W * )U = -iiH(xT + H + H 1 + 5!&))' <14) 

where n d -i = n 2 , n d _i), q 2 = n\ + n\ + ... + n d _ 1 and a 2 = ( — ^p 2 - ^j. Note that in Eq. 
ffT4|) we are using that Cd(s)| s =o — 0. Using the following identity [Bi 



we can see that the first term in the right hand side of Eq. (1141) gives a divergent contribution. 
To proceed we use another useful identity given by 

00 1 1 / 2 \ 

JO, e^ + (2nn d y = 2e{ 1 + ^i)- (16) 

Using both identities given by Eq. f|T5|) and Eq. f|T6|) . it is possible to express the double summation 
that appears in Eq. (fl4"|) by a single summation given by 



t t ln((^) 2 + (2^) 2 )=2 t / ^Q + ^+ai, (17) 

where «i = J2 < H d _ 1 =i S^=-oo l n (l + (2vr?2d) 2 ^. Carrying out the 6 1 integration, we finally arrive 
that Eq. (|T71) can be written as 

where a 2 = «i — S^_ 1= i f 1 + 2 ln(l — e -1 )^) . Since this divergent contribution a 2 is /3-independent 
we will see that it can be eliminated using the third law of thermodynamics. The first term on 
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the right side of Eq. (1181) is a divergent contribution, corresponding to the zero-point energy term. 
Using the following mathematical result [SS] [HE] given by 



n 1+ 



n 2 + b 2 



sinh 2 (7rv / a 2 + b 2 ) 
sinh 2 (7r b) 



(19) 



we can write the last term of Eq. (|14p in a more manageable way. Using the Eq. (fT8l) and Eq. 
ffl9|) . the derivative of the spectral zeta-function in s = can be rewritten as 



d 
ds 



Cd(s) 



s=0 



-2 E 

"<j-i=i 



In 



/sinh(g^+g)\ , , ^ 



— «2 



(20) 



It is possible to show (see appendix A) that, in the finite temperature case, the independent-value 
generating function Qo(<J, h) satisfies Qo(a, h)\ h=g=0 = 1, and 



d 2 
dh 2 



Qo{&, h)\ h -_ 



:<T=0 



r(f) 



2pgS(p\) 



(21) 



In the next section we show that it is possible to obtain a quantum bound in the spatially bounded 
system defined by a self-interacting scalar field in the strong-coupling regime in high temperatures. 
As we will see, for the cases of intermediate or low temperatures, the sign of the renormalized 
zero-point energy is crucial for the validity of a quantum bound for the specific entropy. 



4 The specific entropy for strongly coupled (go ^p v )d theory 

In this section we compute the specific entropy J; of the system. For simplicity, let us define 
In Z(/3, Q, h)\ h=o — lnZ(/3, f2). From Eq. (jBJ) and Eq. ([7j), and using for simplicity that the mean 
energy E(/3, Q) = E and the entropy S(fl, Q) = S, the specific entropy is given by 



S 



— = p-hiz(p,n)[—hiz(p,n) 



E 



d 



dp 



Substituting Eq. ([20) and Eq. (ED) in Eq. (JED we have that In Z(0,£l) is given by 



r(f) 



lnZ{P,tt) = - 
where a' = a + a.i and the quantity h(P) is given by 



\ + UP) 



UP)= E 



In 



/sinh(ff^T 



+ In 1 - e 



T L £g\ ii P q 



+ 



2L 



(22) 



(23) 



(24) 



Defining C\ and C 2 = — that depend only of p and go and do not depend on (3 as 

d = , (25) 

4p(p!)5^o 

the quantity In Jl) can be written in a general form as 

\nZ((3,n) = C 1 -C 2 I 2 ((3). (26) 

It is worth to mention that the quantity C\ corresponds to a divergent expression, C 2 is finite and 
the summation term in the right-hand side of Eq. (1201) is proportional to the zero-point energy. In 
order to renormalize In Z(f3, Q) we first can use the third law of thermodynamics. The derivative 
of lnZ(/9,f2) with respect of (3 yields 

^lnZC9,n) = -C7 2 Aj 2 ( /3 ) j (27) 
where the derivative of I 2 (f3) with respect to (3 is given by 



^ n d _i=X \ \ / e ^ ± 

Using the definition of the mean energy given by Eq. ([6]), Eq. ( [271) and Eq. (1251) we have that the 
unrenormalized mean energy is given by 

£J(v^- h (|f v ^)_ 9COth (^)) + - 3 |^ +9 ) . (29) 

For the case a = 0, i.e., the massless case, using the Eq. (1231) we get that the unrenormalized 
mean energy E(/3,Q) becomes 

7iT(*) °° / 2o \ 
4p(p!)5 c/ p L «d_i=i Ve i - 1 / 

The formula above has the simple interpretation of being phase space sums over the mean energy 
of each mode, where the zero-point energy is included. Note that the divergence that appear in 
the mean energy given by the Eq. (1301) is coming from the zero-point energy, which is given by 

00 

£0 = ^ E (nl+nl + .-.+nU) 1 *, (31) 

f*d_i=i 
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and its sign is given by the ratio between the first and the second terms of the right-hand side 
of Eq. (I3"0"j) . for a negative zero-point energy. The sign of the renormalized mean energy will be 
discussed in details later in this section. 

Let us briefly discuss how to find the renormalized zero-point energy. An analytic regularization 
gives the renormalized zero-point energy. Using the definition of the Epstein-zeta function given 
by 

oo 

A(ai,aa, ...,a k ;2s) = (^nj + a 2 n 2 2 + ... + a k n 2 k )~ s , (32) 

we can find the analytic extension of the Epstein-zeta function in the complex plane, in particular 
at s — — to define the Casimir energy. The structure of the divergences of the analytic extension 
of the Epstein-zeta function is well known in the literature [10] [67] [68] [69]. The renormalized 
zero-point energy E£ is defined as the finite part of a meromorphic function that possesses simple 
poles. Details of this calculations can be found in refs. [H] [70]. As stressed in the Ref. [63j there 
is an ambiguity in the renormalization procedure to find the Casimir energy. 

It is possible to prove that, in a regularization of an ill defined quantity, if the introduction of 
a exponential cut-off yields a analytical function with a pole in the origin, the analytic extension 
of the generalized zeta-function, or the zeta-function method, is equivalent to the application of 
a cut-off with the subtraction of the singular part at the origin [71] [72]. Once we accept the 
advantage of the zeta-function method over the cut-off method with the subtraction of the polar 
terms, due to the fact that it is an analytical extension method, we face a problem: the non- 
trivial scaling behavior of the vacuum energy. This follows from the fact that in an analytic 
extension method it is necessary to introduce a mass parameter /i, i.e., a normalization scale to 
keep the Epstein-zeta function dimensionless for all values of s. If we consider a change in the 
normalization scale \i — > //, it is possible to show that in a (^-dimensional spacetime, the ambiguity 
of the renormalized zero-point energy is related to the B± coefficient by the expression: 

2 

W) = ^oM-^^ln(^). (33) 

Although, in general situations, there is an ambiguity in the renormalization procedure, in our case 
(d{s) \ s =o — Bd = 0. Therefore there is no scaling in the theory and consequently the renormalized 

zero-point energy does not depend on the renormalized scale fi. 

_ 2 

Note that although in the expression for the renormalized mean energy, up to the order (go) p , 
the coupling constant appears, we are interested only in the ratio 4, and the dependence of the 
coupling constant disappears. It is important to stress that this happens only because in the 

_ 2 

strong-coupling expansion, up to the order (g ) p , we are able to split In Z into two contributions: 
the first one proportional to the spectral-zeta function and the second one that has a contribution 
from the independent- value generating function Qo((T, h). 

After this discussions we are able to present the entropy of the system. Substituting Eq. ( 1261) 
and Eq. ( 1271) in the definition of the entropy given by Eq. (JTJ) , we have that the entropy of the 
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system can be written as 

8 = Cl -,K,(!f-±M). (34) 

A system with a unique ground state corresponds to a state of vanishing entropy at zero temper- 
ature. For systems with degenerate ground states this property is also valid if the entropy is a 
extensive quantity Since at zero temperature the system goes to a non-degenerate ground state, 
the entropy must go to zero. The expression of the entropy given by Eq. (|34|) must satisfy the third 
law of thermodynamics, i.e., the entropy of a system has a limiting property that lim^oo S = 0. 
To proceed, lets analyze the limit given by 

hm -=^- J - = hm — I 2 (3) = > , 1 > q. (35) 

Substituting Eq. (133|) in Eq. (I3"4"|) . and using the third law of thermodynamics, we get 

Hm S = Ci = . (36) 

Therefore the first step to find a finite result for lnZ(3,Q) was achieved, since we were able to 
renormalize C% to zero using the third law of thermodynamics. After this step we have 

\nZ(f3,Q) = -C 2 I 2 (3). (37) 

Substituting Eq. (1371) in Eq. (|22|) we can see that for the case a = 0, i.e., the massless case, the 
quotient Ji yields 

1 = 27^(0, (38) 

where we are defining the dimensionless variable £ given by £ = 3/L. Since the field is confined in 
a hypercube, the radius of the smallest (d — l)-dimensional sphere that circumscribes this system 

should be given by R = \ J{d — 1) L. The function T^(£) defined in Eq. ( 1381) is given by 

M4) = — n — t — m — „ ' ^ 



where ej^ = LEq and the positive functions Pd(Q and Rd(Q are defined respectively by 

«d-i=i 

and 



P d (0= E 7rg(e*««-l)~ (40) 



ifc(£) = - E ln(l-e-^) . (41) 



«d-i=i 
12 



Now let us study the function T d (C) given by Eq. (1551) . The quantum bound holds whenever 
Td{0 < 1 for all values of £. From the definition of the function T^(£), given by Eq. (1391) . we have 
that Xrf(£) has a divergent value only if the renormalized zero-point energy is negative. For the 

point £ = £o which satisfies ej + Pd{£,o) = 0, the quantum bound is invalidated. 

Numerical calculations can help us understand the quantum bound. In the Fig. (1) we present 
the plot of the function T^(^) in the case of d = 3 over the interval < £ < 2. Since the 
renormalized zero-point energy is positive [70] . the function T^(£) is also positive for all values of 
£. There is a maximum for some value of £ that we are calling C, ma x, which is near one. For this 
case there is a quantum bound. In Fig. (2) we present the function in the case of d = 4 over 

the interval < £ < 2. Since in this case the renormalized zero-point energy is negative, we have 
that for some value of £ = £o, the function T^(^) diverges. There exists a critical value £ c where, 
for £ > £ c , the specific entropy is unbounded above. 

Let us analyze two cases. The first one is when the renormalized zero-point energy is positive 
(see Fig. 1) and a maximum value for T^(£) appears. The second case, with a negative renor- 
malized zero-point energy, invalidate the quantum bound. For even dimensional spacetime, the 
renormalized zero-point energy is always negative. For the odd dimensional case, it is known that 
for d < 29 this quantity is positive and for d > 29 it changes the sign [TT] . 

For the cases of positive renormalized zero-point energy, an equation for the maximum value of 
Td{0 can be found. The equation for the maximum is given by Rd{£,max) = Q £,max- Substituting 
this i m ax in Eq. (139|) we can find that T d (t; m ax) = J^^T ' Using the same procedure in Eq. (l38j) 
we get -| = P m ax, where f3 max = L£ max . Therefore we can conclude that for odd space-time 
dimensions d < 29 there exists a maximum value for the function Td(C,). 

We can see that the maximum value of Td(C.) depends on the renormalized zero-point energy, 
where for the case d = 3 is less than one. To prove that, for odd d < 29, we have that T^(^) satisfies 
the inequality 7d(£) < 1, let us define an auxiliary function R d (0 that satisfies Rd{£,) < R'di.0- 
This function is given by 

R' d (0 = -[ dQd-i r drr d ~ 2 ln(l-e-^ r ) , (42) 
Jn R Jo v J 

where the angular domain of integration correspond to the region where > 0. Performing 
this integral [66] we have that 

R' d (0 = S d -! T(d - 1) ((d) . (43) 

where the angular term is S d -i = 2 }^^ d -iy Using the Eq. (14"51) in the equation for the maximum, 
i.e., Rdi^max) = &d ^ £,max, we can find that ^ max < £,' max , where 

e ( 2 ^-^ V (44) 
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and we have that T^ max ) < 7r ^^ 1 • In the table 1 we present the maximum values for d = 3 until 
d = 29 for odd d's. 



d 


3 


5 


7 


9 


11 


13 


>) 
£ d 


4.1 x 1(T 2 


6.2 x 10~ 3 


1.1 x 10~ 3 


2.2 x 10~ 4 


4.4 x 10~ 5 


9.4 x 10~ 6 


Td(jiiriax) ^ 


0.3763 


0.2645 


0.2303 


0.2130 


0.2025 


0.1953 




d 


15 


17 


19 


21 


23 


25 


>) 
£ d 


2.0 x 10~ 6 


4.5 x 10~ 7 


1.0 x 10~ 8 


2.2 x 10~ 8 


5.0 x 10~ 9 


1.1 x 10~ 9 




0.1901 


0.1861 


0.1829 


0.1804 


0.1784 


0.1769 



d 


27 


29 


31 


Jr) 
E d 


2.3 x 10~ 10 


3.0 x 10~ n 


-1.1 x lO" 11 


Tdijimax) ^ 


0.1761 


0.1781 


no maximum 



.175 
.15 
. 125 

T d (« 0,1 

. 075 
. 05 
. 025 
0. 



Figure 1: Td(£) as a function of £ for the case of positive renormalized zero-point energy for d = 3. 

Until now we studied the quantum bound for general dimensions based on the summations 
given by Eq. fj4*0l) and Eq. (jUJ). Nevertheless we can find an upper bound function T' d (£) of the 
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0.5 



1 . 5 



T d (£) 



; — -J 











0.5 



1 . 5 



Figure 2: T d (£) as a function of £ for the case of negative renormalized zero-point energy for d = 4. 



function 7d(£) which is more manageable. For this purpose, in a similar way as we have defined 
the function R' d (£), let us define also the auxiliary functions P d (£) and P d (£), that satisfy 

P d (0<K(0, (45) 

and 

p d (o>pm, (46) 

so that the specific entropy satisfies the following inequality 

| < 27^(0, (47) 

where 

rpi ( t\ 1 Z, p 'd(0 + P'd(0 i A o\ 

1 d\s) = — a — T — (r) • ^ 

Without loss of generality we can choose as the auxiliary functions P' d (£) and P d (£) by the 
integrals 

P' d (0 = vr / dQ d ^ r drr*- 1 (e^ r - if 1 (49) 
Jn R Jo v ' 
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and 



P' d (0 = vr cM^J™ drr d - 1 (e^ r - l) 1 . (50) 
Performing these integrals [66], we obtain that P'JJi) and P'^iO are given by 

P' d (0 = irS d ^T(d)ad)(^-) d (51) 



and 



P'>(0 = nS d ^(T(d)ad) - /(d)) (-!-)" , (52) 



where the series f(d) is given by 



B, 



To obtain an upper bound for the specific entropy in a generic Euclidean d- dimensional spacetime 
we have only to substitute Eq. (1151) . Eq. Q5I]) and Eq. ([52]) into Eq. (JUJ). We have that 



fei(d) 

where 



r' d (0 = (or: /, . (54) 



*n(d) = _%L=c(d) ( r(d) + r(d - 1) ) , (55) 



T[ d \fd 
and 

fca(cO = ^(r(d)C(d) -/(«*))• (56) 

It is interesting to study the behavior of the specific entropy for low and high temperatures. 
For the case of high temperatures, we get 

S hAd) , . 

* <M W' (57) 

At high temperatures the dimension in the imaginary direction shrinks to zero and the system 
behaves like a classical system in (d— 1) dimensions where quantum fluctuations are absent. This 
behavior of the specific entropy increasing with (3 in the high-temperature limit was obtained 
by Deutsch in Ref. [21 J . Bekenstein using the condition (3 <C R (high temperature limit) also 
obtained the same behavior in Ref. [18]. Since the thermal energy can compensate the negative 
renormalized zero-point energy, the quantum bound holds. 
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When considering the low temperature behavior of the specific entropy, we can see that the 
problem of the sign of the renormalized zero-point energy can invalidate the quantum bound. In 
this limit we have 



Although some authors claim that the energy of the boundaries of such systems can compensate 
the negative renormalized-zero point energy yielding a net positive energy, for us this is still an 
open question that deserves further investigation. Note that although our results are based in 
a quite particular choice of the shape of the macroscopic boundaries that confine the field in 
the volume Q, it is tempting to think that, using some results from spectral geometry [73], the 
quantum bound could be generalized to more general geometries. 

5 Conclusions and perspectives 

In this paper we studied self-interacting scalar fields in the strong-coupling regime in equilibrium 
with a thermal bath, also in the presence of macroscopic boundaries. In the strong-coupling 
perturbative expansion we may split the problem of defining the generating functional in two 
parts: how to define precisely the independent-value generating functional and how to go beyond 
the independent- value approximation, taking into account the perturbation part. The presence of 
the spectral zeta-function allow us to introduce the boundary conditions in the problem. Using the 

_ 2 

Klauder representation for the independent- value generating functional, and up to the order (g ) * , 
we show that it is possible to obtain a quantum bound in the system defined by a self-interacting 
scalar field in the strong-coupling regime. We established a bound on information storage capacity 
of the strong-coupled system in a framework independent of gravitational physics. 

We have shown that, in the strong-coupling regime, at low and intermediate temperatures 
(/3 ~ L), the quantum bound depends on the sign of the renormalized zero-point energy given by 

E^ r) . For even spacetime dimensions d and also for odd values satisfying the inequality d > 29, 

(r) 

Eq is a negative quantity. Therefore the quantum bound is invalidated. For odd values of 
d, satisfying the inequality d < 29, Eq is a positive quantity. In this situation the specific 

entropy satisfies a quantum bound. Defining as the renormalized zero-point energy for the 
free theory per unit length, we get the following functional dependencies. For low temperatures we 
get Jj < 2irR ( ^^ d J_ 1 , where R is the radius of the smallest sphere circumscribing the system. For 

the case of high temperature, we get that the specific entropy always satisfies a quantum bound, 



Before finishing, we would like to discuss whether the additive normalization energy is accept- 
able for an entropy-energy bound, in respect that at principle we expect that such rate must be 
independent of any additive normalization. Let us discuss briefly the normalization condition, 
since we have the ambiguity in the finite part of the renormalized zero-point energy. As we dis- 




l-d 




given ly|<Ma{. 
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cussed, a merit of the zeta function method over the cut-off method is the fact that it is an analytic 
extension method, where we should introduce a mass parameter to keep the generalized zeta func- 
tion of the problem a dimensionless quantity for all values of s. We would like to stress that this 
is a general feature of any method which is based in the principle of analytic continuation. The 
conclusions of the argument become obvious. We have to consider the effect of a change in the 
normalization scale. Therefore we have scaling properties. These scaling properties are given in 
Eq. f)12p . We proved that the spectral zeta-function in s = is zero and from Eq. fll3p we have 
that C( s ) U=o = Bd. The conclusion is that in the compact region that we are confining the field 
B± = 0. Now we are able to discuss the ambiguity in the finite part of the renormalized quantities 

2 

which belongs to the free theories, i.e., the zero-point energy. As we discussed, if we consider 
a change in the renormalization scale fi, it can be shown that the Casimir energy defined in a 
ci-dimensional spacetime associated with the scale \x and the scale // are related by the expression 
given by Eq. (1551) . Since Bd = 0, we can claim that our results are presented in a normalization 

independent way. 

We would like to stress that without a proof that the spectral zeta-function in s = is 
zero, the correction coming from the boundaries and interaction at principle are not present in a 
normalized independent way, since we have an ambiguity in the zero-point energy. This ambiguity 
was discussed by many authors in distinct situations. For example, for the case of massive scalar 
field theory in a classical background field, Bordag, Mohideen and Mostepanenko [8] claim that, 
in the limit of infinite mass, quantum fluctuations must vanish, so the renormalized energy must 
vanish as well. This was also discussed in refs. [71] [75]- In a four dimensional spacetime, in the 
massless case, it was shown that there is no general normalization condition, if the B 2 coefficient 
is non-zero. 

Although in the framework of analytic extension procedures there is no clear resolution at the 
present to the ambiguity for the zero-point energy if Bd ^ 0, the global energy of the force between 

the hyperplanes has an unambiguous value. We should note that using an alternative procedure, 
the ambiguity of the additive normalization can be fixed in the following way: in the regularization 
procedure, using a exponential cut-off, with subtraction of configurations, if the boundaries go to 
infinity, the renormalized zero-point energy must be zero. Important arguments supporting this 
procedure are based in the demonstration that the expected value of the energy-momentum tensor 
in the vacuum state (since it is a state belonging to discrete spectrum, normalized to unit and 
Poincare invariant) should be zero to ensure that the correct commutation relations of the Lie 
algebra are satisfied by the generators of the Poincare group [76]. The conclusion is that, even 
after obtaining a finite result for the vacuum energy, we still have to use a physical argument 
to fix the value of the energy for some configuration. Therefore our results are acceptable for 
an entropy-energy bound, in respect that such rate must be uniquely defined after raising the 
ambiguity in the finite part of the renormalized zero-point energy. 

There are some continuations for this paper. Still using the scalar field, the Bekenstein bound 
should be investigated assuming more general boundary conditions over the macroscopic bound- 
aries that confine the field. Another interesting point is to investigate the Bekenstein bound in 
different quantum field models still using the strong-coupling expansion. For an alternative method 
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to investigate the strong-coupling regime in quantum field theory, see for example Ref . [77j . The 
generalization of the Bekenstein bound in a theory with vector and spinor field is under investi- 
gation by the authors [75] . 



A Appendix: The Klauder representation for the indepen- 
dent value generating functional 

To give meaning to the independent value generating functional Qo(cr, h), we are using the Klauder's 
result as the formal definition of the independent-value generating functional derived for scalar 
fields in a d- dimensional Euclidean space. It is possible to show that the independent-value gen- 
erating function can be written as 



Q (a,h)=exp^~ J d d x jT ^ (1 - cos(hu)) exp (~a m 2 u 2 - ^ uty . (A.l) 



There is no need to go into details of this derivation (see Ref. |Sj). We would like to point out 
that in Klauder's derivation for the free independent-value model a result was obtained which is 
well defined for all functions which are square integrable in R n i.e., h(x) G L 2 (R n ). Since we are 
assuming that h = cte, we have to normalize our expressions. In order to study Q (a,h) let us 
define E(mo, a, go, h) given by 



/°° du ( 1 9o \ 
- — r (1 — cos(W)) exp — a ml u 2 -h p . (A. 2) 

-do \u\ V 2 p\ ) 



Using a series representation for cosx and using the fact that the series obtained (J2T=i c k fk( u )) 
not only converges on the interval [0, oo), but also converges uniformly there, the series can be 
integrated term by term. It is not difficult to show that 

E(m , a, g , h) = 2 £ ( -^-h 2k f°° du u 2k ~ l exp (-±crm 2 u 2 -^ uA . (A.3) 

Now let us use the fact that the a parameter can be choose in such a way that the calculations 
becomes tractable. Let us choose a — 0. Therefore we have 

E(m ,a,g ,h)\ a=0 = 2j2 [ -7^-h 2k duu 2 ^ 1 exp(-^ u"). (A.4) 

, i 2kl Jo v\ 



k=l 



Let us use the following integral representation for the Gamma function [65] 



f 

Jo 



rf« v - 1 exp(-/ii p ) = -fTp T (- ) , Re(/i)>0 Re(u) > p > 0. (A.5) 
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At this point it is clear that the (go (p p ) theory, for even p > 4, can easily handle applying our 
method. Using the result given by Eq. (1A.5[) in Eq. (IA.40 we have 

00 h 2k 

E(m ,a,g ,h)\ (T=0 = J29(P,k)-2k, (A.6) 



where the coefficients g(p, k) are given by 



2(-l) k 2k 2k 
9(P,k) = - [ -^(p\)TT(-). (A.7) 
p [2k)\ p 



Substituting the Eq. ( 1A.6I) and Eq. ( 1A.7I) in Eq. ( lA.lj) we obtain that the independent-value 
generating function Qo(o", ^)| a=o can be written as 



Qo(a, h)\ a =o = exp 



I f/3 r 00 tfk 

2np Jo j k=1 



(A.i 



It is easy to calculate the second derivative for the independent- value generating function with 
respect to h. Note that Qo(cr, h) \ h=c=0 = 1. Thus we have 

d 2 ( 1 00 h 2k ~ 2 \ ( 1 00 h 2k \ 

Q (a,h)\ a=0 = I -- X] fl'fe k)(2k)(2k - 1) — I exp I — - ^ g(p, k)—^- J +G(g ,p,h), 

(A.9) 

where G(go,p, h) is given by 

(00 l2/c+2<j-2 \ / 1 00 /?, 2fc \ 

£ a(p» fc > g)— 20^- ex p ~o Z g(p» L ( A - 10 ) 

and g(k,q) = k q g(k)g(q). We are interested in the case h = 0, therefore the double series does 
not contribute to the Eq. flA.9j) . since Umh^oG(h) = 0. Using the fact that we are interested in 
the case h = 0, we have the simple result that in the Eq. ( 1A.9I) only the term k = 1 contributes. 
We get 

3 2 r(2) 

^ Qofo &) I h=a=o = r^— ■ (A.H) 

^ 2p<rf(p!)5 

B Appendix: Proof that the value of the spectral zeta- 
function in the origin vanishes, i.e., Cd(s)| s =o = 

As we discussed before, to take into account the scaling properties we should have to introduce an 
arbitrary parameter \x with dimension of a mass to define all the dimensionless physical quantities 



20 



and in particular make the change 

1 d _ , . , 1 d _ , . , 1 / 1 \ . . . 

5S &MU - 5S 0»WU - j m (j^JCdWU. (B.i) 

In this appendix we have a proof that the spectral zeta-function in s = is zero, consequently 
there is no scaling in the theory. The Epstein zeta-function is defined by 

oo , ^ 

Z p (ai,...,Op;2s)= ((ai n x ) 2 + ... + (a p n p ) 2 ) , (B.2) 

m,...,n p =— oo 

where the prime indicates that the term for which all n» = is to be omitted. This summation 
is convergent only for 2s > p. Nevertheless, we can find an integral representation which gives 
an analytic continuation for the Epstein zeta-function except for a pole at p = 2s [ID] . This 
representation is given by 

(vr?7) _s r(s) Z p (ax,...,a p ;2s) = 

— + — ^7- + rT s / dxx"~ l (0(0, ...,0;a 2 £, ...,a 2 x) - l) 
s p — 2s J n v p ' 

+ri Q*-p)/2 J™ dxxtr-W 2 - 1 (f(0, ...,0;x/a 2 , x/aj) - l) , (B.3) 

where 7/ P//2 is the product of the p's parameters a; given by r/ p / 2 = a\...a p , and the generalized 
Jacobi function ~&(zi, ...,z p ;xi, ...,x p ), is defined by 



&(zi, ...,z p ] xi, x p ) = JJ d(zi\ Xi) , (B.4) 

i=i 

with d(z; x) being the Jacobi function, i.e., 

oo 

a;) = ^ e" (2n2 ~ n2;c) . (B.5) 



n=— oo 



Using this integral expression for the Epstein zeta-function, given by Eq. (1B.3|) . we can find that 

Z p (ax,...,a p ;2s)\ s=0 = -1 , (B.6) 
for any p > 1. To proceed, let us define the function Z^(ai, ...,a p ; 2s), given by 

oo oo _ 

Zto(a 1 ,...,a p ',2s)= £ £ ((oi m) 2 + ... + (a p n p ) 2 )" . (B.7) 

»i,...,n 9 =l n q +i,...,n p =— oo 
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Using the result given in Eq. (1B.6[) we can show that, after performing the analytic continuation 
of the function Z^ q '(ai, ...,a p ; 2s), the following property holds 

ZW(ai,...,o p ;2 S )|. =0 = 0, (B.8) 

where this result is valid only for < q < p. We can prove this property by induction. First, let 
us verify that for = 1 the above property hold. Therefore, assuming that is valid for q, we have 
only to show that is true for q + 1. For q = 1 we have that 

Z p (ax, a p ; 2s) | s=0 = Z p (a 2 , ...,a p ; 2s)| s=0 + 2Z ( p l \a 1 , a p ; 2s)| s=0 . (B.9) 

Since p > 1 we can use the property given by Eq. flB.6j) . for the two first terms of Eq. (IB .91) 
and verify that Z^(ai, a p ; 2s) | s =o — 0. The next step in the proof by induction is to assume 

the validity of this property for some q, i.e., Z^ q '(ai, a p ; 2s) | s= o = 0, with p being arbitrary, 
but satisfying the condition < q < p, then we must verify the validity of this property for 
q + 1, i.e., Zp, +1 \ai, a p i ; 2s) | s=0 = with p' also being arbitrary but satisfying the condition 
< q + 1 < p 1 . From the following property 

Z^){ai, ...,Op> ; 2s)| s=0 = Z ( v q )_ x (a u ...,a q ,a q+2 , ...,a p >; 2s)| s=0 + 2Z^ +1) (ai, a p /; 2s)| s=0 , (B.10) 

since < q < p' — 1 and using the assumption of the validity of this property for arbitrary q, 
given by Eq. ( IB. 81) . we can see that the two first terms in Eq. (IB. 101) vanish. Therefore we finally 

proved that Z^,(ai, a p > ; 2s) | s=0 = 0. We are interested in a particular case of this property, 
given by 



Zf 1) (ai,...,a P ;2s)| a=0 = ( ((^i ^i) 2 + ••• + (« P n p) 2 

Lni,...,n p _i=l n p =-oo 



= 0. (B.ll) 

s=0 
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